Magnetic ordering and magnetodielectric phenomena in CoSe04 
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CoSe04 has a structure consisting of edge-sharing chains of Co^^ octahedra which are held to- 
gether by Se04~ tetrahedra via shared oxygen atoms at the edges of the octahedra. DC magnetiza- 
tion measurements indicate a transition to an ordered state below 30 K. Powder neutron diffraction 
refinements suggest an ordered state with two unique antiferrromagnetic chains within the unit cell. 
Isothermal magnetization measurements indicate a temperature-dependent field-induced magnetic 
transition below the ordering temperature. From neutron diffraction, we find this corresponds to 
a realignment of spins from the canted configuration towards the c-axis. The dielectric constant 
shows a change in slope at the magnetic ordering temperature as well as a quadratric dependence 
on the external magnetic field. 

PACS numbers: 75.50.Ee,75.30.Kz 



I. INTRODUCTION 

Magnetodielectrics are materials in which the dielec- 
tric properties couple to changes in the magnetic or- 
der. The ability to find and design new single-phase 
materials which exhibit this coupling between the spin 
and charge degrees of freedom has significant techno- 
logical implications in the development of magnetic sen- 
sors and field- tunable dielectrics.^^ Systems showing in- 
commensurate and noncollinear spin structures such as 
CoCr204^'^, MnsOj^, SeCuO^l, and NisVsOg'^ have all 
been extensively studied to understand the nature of 
such interactions. Recent work characterizing the mag- 
netoelectric coupling in the magnetic chain compounds 
MnW04!^and LiCu2 02-'^ indicate that systems with 
reduced crystallographic dimensionality should be of par- 
ticular interest because of the strong interplay between 
the charge, lattice, and magnetic degrees of freedom that 
occurs in these materials) Reduced dimensionality 
has been found to result in noncollinear and canted anti- 
ferromagnetic spin structures which break symmetry re- 
strictions and allow electric polarization to develop. 

Here we discuss the preparation and detailed mag- 
netic and dielectric characterization of CoSe04; a mag- 
netic chain compound which adopts a crystal structure 
analagous to /3-C0SO4 (Fig. [T]). The orthorhombic unit 
cell contains chains of octahedrally coordinated Co^+ 
which are bound together by Se04~ tetrahedra via shared 
oxygen atoms at the edges of the octahedra. The chains 



of octahedra are tilted with respect to neighboring chains 
with the axis of the octahedra tilted —35° and +35° off 
the a-axis for the chains on the edge and in the center 
of the cell respectively. Previously reports have charac- 
terized the nuclear and magnetic structure using neutron 
diffraction,' ^^ ' -^^ but no reports on the magnetic suscep- 
tibility or dielectric properties have been made to date, 
and the effects of magnetic fields have not been fully in- 
vestigated. Using detailed magnetization, heat capacity, 
powder neutron diffraction, and dielectric measurements, 
we show that the transition to long-range antiferromag- 
netic below 30 K is accompanied by a signficant change 
in the slope of the temperature dependence of the di- 
electric constant. We also find a field-induced change 
in the magnetic structure. The dependence of the di- 
electric constant on the external field, which is quadratic 
in nature, changes dramatically above the field-induced 
magnetic transition. 



II. EXPERIMENTAL METHODS 

The title compound was prepared following a previ- 
ously reported procedure^^ by neutralizing a solution of 
selenic acid (2.559 g 23.06 mmol) of Se02 dissolved in 
50 cm^ of water) with Co2C03(OH)2 (1.629 g, 7.69 mmol) 
at 70° C and recrystallizing in water. The final crystal- 
lization step was performed by allowing the remaining 
water to evaportate at room temperature over several 
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FIG. 1: (Color Online) Illustration of the crystal structure of 
CoSe04. Chains of octahedrally coordinated Co^^ (purple) 
are bound together by Se04 tetrahedra via shared oxygen 
atoms at the edges of the octahedra. Note that the chains 
are alternately tilt +35° and —35° off the a-axis. (a) View 
down the c-axis of the unit cell, (b) View down the c-axes to 
illustrate the chains of edge sharing CoOe 



days. Faceted crystals of dark red CoSe04 -61120 up 
to 2 cm on a side were obtained from this proceedure. 
These crystals were ground and dehydrated at 125° C 
overnight to yield a bright pink powder of CoSe04-H20 
which was subsequently heated at 315° C for several days 
to produce a light violet powder of anhydrous CoSe04. 
Co2C03(OH)2 was prepared following a separate previ- 
ously reported procedure. Stoichiometric amounts of 
CoN03-5H20 and Na2C03 were separately dissolved in 
H2O with a 2:1 (NaCOsiCoNOa) volume excess. The so- 
lution of C0NO3 was added to the solution of Na2C03 
which had been preheated to 90° C. After stirring for 2 
hours the resulting dark purple precipitate was collected 
by vacuum filtration and subsequently washed with H2O 
and ethanol. Purity of the precursors and final prod- 
uct was confirmed by powder X-ray diffraction on a 
Philips XPERT MPD diffractometer operated at 45 kV 
and 40 mA. 

Temperature dependence of the DC magnetization was 
measured on well-ground powder samples using a Quan- 
tum Design MPMS 5XL SQUID magnetometer. The 
specific heat data were collected using the semiadiabatic 
technique as implemented in a Quantum Design Physi- 
cal Property Measurement System (PPMS), under zero 
applied field as well as under a 50kOe field. The mea- 
surement was made by mixing the compound with equal 
parts by mass of Ag powder and pressing into a pellet 
in order to improve thermal coupling to the stage. The 
contribution from Ag was measured separately and sub- 
tracted. Variable temperature neutron diffraction data 
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FIG. 2: Field-cooled (FC, open symbols) and zero-field cooled 
(ZFC, closed symbols) magnetic susceptibility of a powder 
sample of CoSe04 acquired under increasing strengths of the 
external magnetic field. At 30 K, the system orders to an 
antiferromagnetic state. A weak ferromagnetic component 
is found in the 100 Oe data, suggested from the separating 
ZFC and FC traces. Note also that the cusp which indicates 
the magnetic ordering is smeared out with increasing field 
strength. 



under a magnetic field were collected on the BTl pow- 
der diffractometer at the National Institute of Standards 
and Technology, Gaithersburg, MD using a wavelength 
of 2.08 A. Dielectric properties were measured by at- 
taching polished copper electrodes to opposite faces of a 
cold-pressed pellet of CoSe04, using a Quantum Design 
PPMS for temperature and field control, and an Agilent 
4980A LCR meter to measure the capacitance. 



III. COMPUTATIONAL METHODS 

Density functional theory (DFT) calculations were 
performed using the Vienna Ab-initio Simulation Pack- 
age (VASP)P^^^ at the experimental lattice parame- 
ters. The projector augmented- wave (PAW) methocPl 
was used together with the local density approxima- 
tion (LDA) Ceperley-Alder exchange correlation func- 
tional.^^ A plane- wave energy cutoff of 500 eV and a 
6x3x9 F-centered Monkhorst-PaclJ22 /c-point mesh cor- 
responding to 162 irreducible /c-points was used to sam- 
ple the Brillouin zone. The tetrahedron method with 
the Blochl correction^^ was used for Brillioun zone in- 
tegrations. Correlation was treated using the LDA+/7 
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FIG. 3: (Color Online) (a) Isothermal magnetization loops at 
various temperatures below the Neel temperature. Note the 
field induced transition which pushes to higher fields at low 
temperatures. 



formalism within the rotationally invariant approach of 
Liechtenstein et al?'^ The applicability of the LDA+/7 
formalism is validated in the example of (f Co^+ be- 
cause crystal symmetry results in an absence of orbital 
degeneracy. Several values of U were tested with only 
qualitative changes to the band structure being found. 
Therefore a /7 of 5 eV and a J of 1 eV was chosen to 
agree with previously published DFT results calculated 
on octahedrally coordinated Co^+.^^ 



IV. RESULTS AND DISCUSSION 

The high temperature region (200 K to 300 K) of the 
inverse susceptibility was fit to the Curie- Weiss equation, 
C /{T — Qcw)' A Curie- Weiss temperature of — 36K and 
an effective moment jieff of 4.38 /i^ was determined for 
the data collected under a lOkOe field. The effective 
moment lies between the spin-only value of 3.87 /i^ and 
the value of 5.2 /i^ expected for octahedrally coordinated 
Co2+ {(f, tlgel, S=3/2, L=3) with a fully unquenched 
orbital con tribution as obtained using the relationship 
jjLL+s = V'4S(S + 1)+L(L + 1).1^ The existence of an 
orbital contribution is common from high-spin (f sys- 
tems in an octahedral crystal field due to the orbital 
degeneracy in the t2g levels. The small reduction of the 
orbital contribution likely arises from the irregular octa- 
hedral environment of the magnetic sites reducing some 
hybridization with the surrounding oxygen anions. 

Fig. |2] shows the temperature dependence of the mag- 
netic susceptibility of CoSe04 under a variety of mag- 
netic field strengths. The spins begin to order below 
30 K which is very near the expected Curie- Weiss tem- 
perature indicating that there is very little frustration 
from competing exchange interactions. When measured 
in a 100 Oe field, a sudden jump in the susceptibility 





20 




15 


o 




E 


10 








5 









3 




2.5 




2 


mol" 


1.5 




1 




C/D 


< 


0.5 



i 1 1 1 1 


1 ' 1 ' 1 ' 1 ^ 


" (a) 
- 


n — 




^^jj^'^ o kOe 




° 50 kOe 


1 1 1 1 1 


1 , 1 , 1 , 1 T 




FIG. 4: (a) Temperature dependence of the specific heat of 
CoSe04 measured on a powdered sample in if = OkOe and 
50 kOe field, (b) Entropy released due to the magnetic order- 
ing, obtained by integration of the magnetic heat capacity. 



occurs at 30 K which can be attributed to some weak fer- 
romagnetism which results from a canted antiferromag- 
netic arrangement of the spins. Measuring the sample in 
larger fields minimizes the contribution from the weak- 
ferromagnetic component and gives rise to the cusp-like 
behavior expected from a well-compensated antiferro- 
magnet. This cusp is found to broaden with increasing 
field strength indicating that the magnetic order may be 
driven to a different ground state in the presence of suf- 
ficiently large external fields. To further examine the 
nature of the field-dependence of the magnetic order, 
isothermal magnetization loops were collected at a va- 
riety of temperatures (Fig.[3|. An upturn from the linear 
field-dependence expected for a well-behaved antiferro- 
magnet is found beginning in a field of approximately 
20kOe at 28 K. The upturn shows a strong temperature 
dependence, with the field required to induce it increas- 
ing to 45 kOe by 18 K and subsequently out of the 50 kOe 
range of the magnetometer used here by 10 K. Such a 
field-induced magnetic transition may be attributed to 
overcoming of the magnetocrystalline anisotropy of the 
octahedral Co^+ cations which would manifest as a spin- 
flop transition where the spins realign from one preferred 
axis to another. 

Fig. [4] (a) shows the specific heat of CoSe04 measured 
in zero field, and under a an external magnetic field of 
50kOe. The specific heat shows a sharp anomaly at 
the magnetic ordering temperature of 30 K in the ab- 
sence of an external magnetic field. Application of a 
50kOe field results in a smearing out of the transition 
in a manner similar to the behavior observed in the 
susceptibility data shown in Fig. [2] The lattice contribu- 
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FIG. 5: (Color Online) Densities of state of CoSe04 from a 
collinear spin structure calculation. The magnetic structure 
used in the calculation was that of antiferromagnetic chains 
aligned antiparallel to neighboring chains. 



tion to the specific heat was approximated by fitting the 
data above the magnetic transition to a polynomial ex- 
pansion (/3T^+7T^+^T'^). Subtracting the lattice con- 
tribution and integrating Cp/T yields a change in en- 
tropy due to magnetic ordering of 2.1Jmol~^K~^ and 
2.5Jmol~^K~^ in zero field and a 50 kOe field respec- 
tively. This value is significantly smaller than the value 
of 11.5 J mol~^ predicted by the Boltzmann equa- 
tion (A^ = i?ln(2S + 1), S = 3/2). This considerable 
diminution of the spin entropy from what is expected may 
indicate that the reduced dimensionality of the system 
prevents complete ordering of the spins, or that possibly 
the order within chains is not well correlated with that 
in neighbording chains. 

To better understand the electronic structure and 
chemical bonding in the title compound, density func- 
tional theory calculations were performed using the 
LDA+/7 formalism. Fig.[5]illustrates the density of states 
(DOS) calculated for the experimental room temperature 
structure of CoSe04. The DOS show a significant degree 
of overlap between the Se-p and 0-p states and between 
Co-d and 0-p states below the Fermi energy. However, 
there is a gap of approximately 250 meV between the Co 
and Se states. It can be seen that a majority of the Se and 
O overlap is related to bonding between Se and O3. The 
03 position forms two bonds within the Se tetrahedron 
which are 1.56 A; short in comparison with the Se-01 
and Se-02 bond lengths which are 1.74 A and 1.78 A re- 
spectively. Such bonding is consistent with the idea that 
the Se04 tetrahedra contains two double and two single 
bonds with surrounding oxygens to satisfy its octet. It 
is interesting to note that the Se-03 bonds correspond 
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FIG. 6: Low angle region of the powder neutron-diffraction 
patterns of CoSe04 (BT1,NIST) obtained at different tem- 
peratures. 



to the shared oxygens at the corners of the CoOe octa- 
hedra within a single chain. This arrangement of double 
bonds can be expected to significantly influence magnetic 
superexchange between the CoOe chains. 

The DOS calculated using the LDA formalism results 
in metallic densities of state, but by accounting for cor- 
relations within LDA+[/, with U = 5eV and J = leV, 
a gap of approximately 3eV opens between the occu- 
pied and unoccupied minority spins. We also note that 
the band gap which is approximately 2 eV corresponds to 
the gap between the highest filled states and unoccupied 
0-p orbitals. The distorted octahedral environment of 
the Co^+ breaks the degeneracy of the unoccupied Cg or- 
bitals corresponding to d^^ and dx2_y2 with unoccupied 
dxz above the Fermi energy. 

The antiferromagnetic state is stabilized by more than 
0.5 eV compared to a ferromagnetic alignment of spins. 
The calculated spin magnetic moment on each Co is 
2.72 (iB which is slightly reduced from the full value of 
3 /i^ expected for a d^ cation in an octahedral coordina- 
tion environment. Such a reduction may be associated 
with the covalent bonding to the surrounding oxygens 
or may be related to the fact that the calculation is re- 
stricted to a collinear spin configuration which is not the 
true ground state observed via neutron diffraction. 

The thermal evolution of the neutron-diffraction pat- 
terns collected from 300 K to 5 K are shown in Fig.[6]with 
Table[T| showing a summary of refined parameters. Three 
magnetic reflections appear below 10 K at 12°, 18°, and 
23° , all of which are consistent with a propagation vector 
k = 0. It was only possible to simultaneously fit all three 
magnetic peaks by using the basis vectors of the irre- 
ducible representation Fi. The resulting magnetic struc- 
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TABLE I: Summary of the results of Rietveld structure re- 
finement of variable temperature neutron diffraction data ob- 
tained in zero magnetic field. 





300 K 


20 K 


5K 


a (A) 


9.0421(3) 


9.0355(4) 


9.0344(4) 


6(A) 


6.7669(3) 


6.7397(3) 


6.7398(3) 


c(A) 


4.8860(2) 


4.8759(2) 


4.8757(2) 


M (mb) 




3.14(8) 


3.50(8) 


y (A^) 


298.964(2) 


296.927(2) 


296.886(2) 


Rnnc. (%) 


8.3 


8.8 


8.4 


i^mag. (%) 




15.7 


14.1 



ture is illustrated in Fig.jSj The spins on each Co align 
antiferromagnetically down the length of the chain and 
with respect to the neighboring chains. The moment of 
each spin has components along all three axes of the unit 
cell as listed in Table [TI| with a total magnetic moment of 
3.59 /iB at 5K. The magnetic structure in the absence of 
a magnetic field shown in Fig. [8] agrees well with the pre- 
viously reported structure of Fuess.^ We note that while 
Fuess reported the structure in Pbnm we have chosen the 
standard setting of Pnma. The resulting spin configu- 
ration (compared to that reported by Fuess) has spins 
oriented largely in the ac plane with the moments tilting 
60° (60°) off the a-axis, 70° (69°) off the 6-axis, and 44° 
(41°) off the c-axis again at 5K. The moment within in- 
dividual chains do not cancel completely with a resulting 
uncompensated moment along the a-axis which is can- 
celed by the neighboring chains. This uncompensated 
moment within the individual chains may be a source 
of the weak ferromagnetism observed in small external 
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FIG. 7: (Color Online) Low angle region of the powder 
neutron-diffraction pattern of CoSe04 obtained at 5 K and 
broken down in to the contribution from the nuclear, mag- 
netic, and a small non-magnetic impurity phase of Se205. 



magnetic fields. 

The effect of applying a magnetic field of 70kOe is 
to first reduce the total magnetic moment to 2.60 /i^ at 
5K. The field also causes the spins to reorient with re- 
spect to the unit cell. At 20 K the moments form angles 
of 66° off the a-axis, 64° off the 6-axis, and 37° off the 
c-axis. From examination of the isothermal magnetiza- 
tion measurements in conjunction with the reorientated 
moments the field-induced magnetic transition appears 
to corresponds to a spin-flop transition on to the c-axis. 

Fig. [9] (a) shows the inverse magnetic susceptibility nor- 
malized by the fits of the high temperature region (be- 
tween 200 K and 300 K) to the Curie- Weiss formula. Nor- 
malizing in this manner emphasizes deviations from the 
Curie- Weiss equation in the form of short-range correla- 
tions between the spins which develop above the ordering 
temperature.'^ Fig.|9|b) shows the temperature depen- 
dence of the dielectric constant with a very clear change 
in slope in the dielectric constant at the magnetic order- 
ing temperature. 

To further characterize the coupHng between the mag- 
netism and dielectric properties, the dielectric constant 
was measured as a function of external magnetic field. 
A quadratic dependence of the dielectric constant with 
small fields was found and is shown in Fig.[To| Above 
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FIG. 8: Proposed magnetic structure of CoSe04 as deter- 
mined from Rietveld refinements of the neutron diffraction 
pattern obtained using a A = 2.08 A at 2K. (a) View down 
the 6-axis of the Pnma structure, (b) View down the c-axis. 
The small light grey atoms are oxygen while the larger and 
darker grey atoms are selenium 
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TABLE II: Components of the magnetic moment on each Co 
in the magnetic structure described by Fi. Note that Co at 
the origin occupies a 4a Wyckoff site. The Co within the unit 
cell are: Col (0, 0, 0), Co2 (|, 0, |), Co3 (0, |, 0), Co4 (|, 



2 ' 2^ 





OkOe 


70kOe 




My 


My Mz 


20 K Col 
Co2 
Co3 
Co4 


1.51 1.05 2.51 
-1.51 -1.05 2.51 
-1.51 1.05 -2.51 

1.51 -1.05 -2.51 


0.86 0.93 1.68 
-0.86 -0.93 1.68 
-0.86 0.93 -1.68 

0.86 -0.93 -1.68 


5K Col 
Co2 
Co3 
Co4 


1.76 1.19 2.90 
-1.76 -1.19 2.90 
-1.76 1.19 -2.90 

1.76 -1.19 -2.90 


1.24 1.14 1.98 
-1.24 -1.14 1.98 
-1.24 1.14 -1.98 

1.24 -1.14 -1.98 



the field-induced magnetic transition, the dielectric con- 
stant no longer exhibits a quadratic increase but rather 
changes abruptly to what appears to be a linear de- 
crease. This could reflect the possibility that the low- 
field magnetic structure allows for a magnetodielectric 
response whereas the high-field magnetic structure does 
not. A similar temperature and field- dependent behavior 
has been reported by Nenert and coworkers^^ for a hy- 
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FIG. 9: (a) Inverse magnetic susceptibility normalized by the 
values extracted from fits to the Curie- Weiss formula in the 
high temperature region. Note that this manner of plotting 
emphasizes short-range correlations which cause deviations 
from ideal Curie- Weiss behavior which is illustrated as the 
dashed line, (b) Temperature dependence of the dielectric 
constant of CoSe04 measured on a pressed pellet of a poly- 
crystalline sample at a frequency of 1 MHz. The dashed line is 
a guide-to-the-eye to emphasize the transition at the magnetic 
ordering temperature. 
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FIG. 10: Change in the dielectric constant measured at 1 MHz 
plotted as a function of the strength of the external magnetic 
field at different temperatures. Note the quadratic depen- 
dence until the field-induced magnetic transition 



brid Cr(II) inorganic-organic material. They also find a 
quadratic change with external magnetic field above the 
magnetic ordering temperature as well as below and at- 
tribute the behavior to the always allowed P^H^ order 
parameter. In that aspect, our results differ in that we 
see no response in the dielectric constant in the param- 
agnetic regime. 



Summary 

We have discussed the preparation and detailed mag- 
netic and dielectric characterization of a magnetic chain 
compound. Magnetization and specific heat measure- 
ments show a transition to a long-range canted an- 
tiferromagnetic state below 30 K and a temperature- 
dependent field-induced magnetic transition. Powder 
neutron diffraction shows the spins lay predominantly in 
the ac-plane with the moments tilting 60° off the a-axis, 
70° off the 6- axis, and 44° off the c-axis again at 5 K. Tem- 
perature and field dependent measurements of the dielec- 
tric constant show an anomaly at the magnetic ordering 
temperature, which, at fixed temperature, is quadratic 
in the external magnetic field. Temperature dependent 
measurements of the dielectric constant show an anomaly 
at the magnetic ordering temperature. We also find that 
the field-induced magnetic transition causes an abrupt 
change in the field dependence of the dielectric constant. 
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